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ABSTRACT. Finite differences of values of the Riemann zeta function at the integers are 
explored. Such quantities, which occur as coefficients in Newton series representations, 
have surfaced in works of Bombieri-Lagarias, Maslanka, Coffey, Baez-Duarte, Voros and 
others. We apply the theory of Norlund-Rice integrals in conjunction with the saddle- 
point method and derive precise asymptotic estimates. The method extends to Dirichlet L- 
functions and our estimates appear to be partly related to earlier investigations surrounding 
Li's criterion for the Riemann hypothesis. 



Introduction 

In recent times, a variety of authors have, for a variety of reasons, been led to consider- 
ing properties of representations of the Riemann zeta function £(s) = l/rt s as a Newton 
interpolation series. Amongst the many possible forms, we single out the one relative to a 
regularized version of Riemann zeta, namely, 



(i) coo-^-r = X>i) n &n 

s — 1 ' 

71=0 

where ( S J is a binomial coefficient: 

_ s(s - 1) • • • (s - n + 1) 



Corollary[T]in Section|5]establishes that the representation ([T]) is valid throughout the com- 
plex plane, its coefficients being determined by a general formula in the calculus of finite 
differences (see [13] OH |20) and Section[T]below): 

(2) K = n(l - 7 " H n -i) -~+£(f) (~l) fc C(fc), 

k=2 ^ ' 

Here, H n = 1+ 1 +• • • — is a harmonic number. Although the terms in the sum defining b n 
become exponentially large (of order close to 2"), the values of the b n turn out to be 
exponentially small, while exhibiting a curious oscillatory behavior. We shall indeed prove 
the estimate (Theorem[T]in SectionlH 



(3) b n = ( ^) V4 e~ 2 ^ cos (2V^ - y ) + (V 1 
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Our first motivation for investigating (Q~|) and §2^ was an attempt by one of us, Linas 
[2003, unpublished; available at http : //linas . org/math/poch-zeta .pd.f], to 
obtain alternative and tractable expressions for the Gauss-Kuzmin-Wirsing operator of con- 
tinued fraction theory In particular, Linas' computations at that time revealed that the b n 
tend rather fast to and exhibit a surprising oscillatory pattern, both facts crying for expla- 
nation. The present paper, essentially elaborated in early 2006, represents the account of 
our joint attempts at understanding what goes on. 

Second, the zeta function has received attention in physics, for its role in regularization 
and renormalization in quantum field theory. Motivated by such connections, Maslanka in- 
troduced in [ 16 1 what amounts to the Newton series representation of a regularized verson 
of C(s), namely, (1 — 2s)£(2s). (Further numerical observations relative to the corre- 
sponding coefficients are presented by this author in ifTTI .) The growth of coefficients in 
Maslanka's expansion has been subsequently investigated by Baez-Duarte [2]. In particu- 
lar, Baez-Duarte's analysis implies that the coefficients decrease to faster than any power 
of 1/n. The methods we develop to derive Equation (f3]) can be easily adapted to yield 
refinements of the estimates of Baez-Duarte and Maslanka. 

A third reason for interest in the representation (Q3 and the companion coefficients (O is 
Li's criterion lTT5l for the Riemann Hypothesis (RH). Let p range over the nontrivial zeros 
of C(s). Li's theorem asserts that RH is true if and only if all members of the sequence 



are nonnegative. Bombieri and Lagarias [4] offer an insightful discussion of Li's criterion. 
Coffey [5], following Bombieri and Lagarias [4, Th. 2], has expressed the A„ as a sum of 
two terms, one of which is an elementary variant b n defined by 



Theorem 2 of fH amounts to the property that the coefficients a„(l,2) decrease to 0. 
As we shall see in Section [6] and right after Equation (l48l i in Section [7] the methods origi- 
nally developed for estimating b n yield precise asymptotic information on a n (l, 2) as well. 
Though the sums we deal with count amongst the far easier ones, our precise asymptotic 
estimates may contribute to bring some clarity in this range of problems. 

In this essay, we approach the problem of asymptotically estimating differences of zeta 
values by means of a combination of two well established techniques. We start from a 
contour integral representation of these differences as defined by (O (for this technique, 
see especially Norlund's treatise l20l and the study [9]), then proceed to estimate the cor- 
responding complex integral by means of the classical saddle-point method of asymptotic 
analysis 0ET1- Our approach parallels a recent paper of Voros [25 1 (motivated by Li's 
criterion), which our results supplement by providing a fairly detailed asymptotic analysis 
of differences of zeta values. 

Section[T]presents the construction of a Newton series for the zeta function, and presents 
generating functions for its coefficients. This is followed, in Section [2] by a brief exami- 
nation of numerical results. Section |3]gives the Norlund integral representation for the co- 
efficients, of which Section|4]provides a careful saddle-point analysis. The convergence of 
the Newton series representation (dJ is then discussed in Section|5] and Section|6]develops 
the corresponding analysis for Dirichlet L-functions. We end with a conclusion, Section|7] 





n-l — 



-( 7 -l + log2). 
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outlining other applications of Norlund integrals in the realm of finite differences and zeta 
functions. 

1. NEWTON SERIES AND ZETA VALUES 

This section defines the Newton series for the Riemann zeta that is to be studied, demon- 
strates some of its basic properties, and gives some generating functions for its coefficients. 
In this paper, a Newton series is defined as 



(4) 



OO 

E 

n=0 



(-l) n c n 



Given a function <f>(s), one may attempt to represent it in some region of the complex 
plane by means of such a series. Since the series <£>(s) terminates at s — 0, 1, 2, . . ., the 
conditions (j>(m) = $(m) at the nonnegative integers imply that the candidate sequence 
{c„} is linearly related to the sequence of values {(j)(m)} by 



0(m) = 5^(-l) n c B 



71=0 



The triangular system can then be inverted to give (by the binomial transform ifTTl p. 192], 
or its Euler transform version [22, p. 43], or by direct elimination) 



(5) 



E 

k=Q 



(-l) fe 0(fc), n = 0,1,2, 



This choice of coefficients for (0]l determines the Newton series associated to tp. The 
coincidence of the function cf> and its associated series $ is, by construction, granted at 
least at all the nonnegative integers. The validity of $(s) = cf>(s) is often found to extend 
to large parts of the complex plane, but this fact requires specific properties beyond the 
mere convergence of the series in 



In the case of the Newton series for £(s) 
the coefficients in the form 



(6) 

where 

(7) so 



so-ns 1+ J2 U (-1) 



C(«) 



fc=2 
_ 1 

~ 2' 



l/(s — 1), the general relation (0 provides 
1 



k- 1 



s - 1 

The harmonic numbers appear a^| 

n\ (-l) k 



si = lim 

s— >i 



7- 



(8) 



E v /, 

k=2 



1 — n + nH„ 



Equations (O, (0, © then entail that the b n , as defined by (O, are indeed the coefficients 
of the Newton series associated to ((s) — l/(s — 1). A proof that the equality $(s) = 



'identity (SJ is easily deduced from the classical partial fraction decomposition (see, e.g., 11 II p. 188]), 

n! _ , m\ (-l) fc 

x(x + 1) ■ ■ ■ (x + n) ~ 2-* 

upon letting x —* — 1. 



kl x + k ' 
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C(s) — l/(s — 1) holds for all complex s is given in Section|5] following the asymptotic 
analysis of the coefficients b n and based on a theorem of Carlson. 

Before engaging in a detailed study of the b n , we note a few simple facts regarding their 
elementary properties. Consider the quantities 



(9) 



k=2 



(-i) fe C(fc), 



which represent the nontrivial sums in the definition (f2]i of b„ and are, up to minor adjust- 
ments, differences of zeta values at the integers. Expanding the zeta function according to 
its definition and exchanging the order of summations in the resulting double sum yields 



(10) 



e>i 



This rather simple sum shows a remarkably complex behavior; elucidating its behavior is 
one of the principal topics of this paper. 

The ordinary generating function for the sequence {5 n } is also of interest. Given the 
classical expansion [l] §6.3] of the logarithmic derivative t/j(z) = V (z) /T(z) of the 
Gamma function, 

^(1 + z ) + 7 = C (2)z - C(3)z 2 + C(4)z 3 - • • • , 

one finds, by the usual generating function translation of the Euler transform 1191 p. 311] 
or by an immediate verification based on the binomial theorem: 



(11) 



5>* n = ?r 



n>2 



1> 



1 



+ 7 



The exponential generating function for the sequence {S n } reflects ( fTOb and is even sim- 
pler: 



(12) 



(- z y 



n>2 



- 1 



(The first equality is easily verified by expanding e 2 and expressing the coefficient of z n / n! 
in the product as a convolution, itself seen to coincide with (|9]).) 



2. Experimental analysis 

Detailed experiments on the coefficients b n conducted by one of us are at the origin of 
the present paper and we briefly discuss these since they illustrate some concrete numerical 
aspects of the sequence (b n ) while being potentially useful for similar problems. As it is 
usual when dealing with finite differences, the alternating binomial sums giving the b n 
involve exponential cancellation since the binomial coefficients get almost as large as 2 n . 
We started by conducting evaluations of the b n up to n ~ 5000, which requires determining 
zeta values up to several thousand digits of precision. Note that the zeta values can be 
computed rapidly to extremely high precision using several efficient algorithms (e.g., O), 
which are available in symbolic computation packages (Maple) and numerical libraries 
(Pari/Gp). 
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A quick inspection of numerical data immediately reveals two features of the con- 
stants b n : they are oscillatory with a slowly increasing (pseudo)period and their absolute 
values are rapidly decreasing. For instancqj 

b l = -7.72156- 1CT 2 , b 2 = -9.49726 • 10~ 3 , b 5 = +7.15059 • 1(T 4 , 



"5 

&io = -2.83697 ■ 10~ 5 , b 2Q = +2.15965 • 10" 9 , 6 5 o = -1.08802 • 10" 11 . 



A numeric fit of the oscillatory behavior of the function was first made. There are sign 
changes in the sequence {&„} at 

n = 3, 7, 13, 21, 29, 40, 52, 65, 80, 97, 115, 135, 157, 180, 

the values growing roughly quadratically. A good fit for the kth sign-change was found 
experimentally to be of the form q(k) — jk 2 + 0(fc). The quadratic polynomial is then 
easily inverted to give an approximate oscillatory behavior of the b n . Once the oscillatory 
behavior had been disposed of, the task of quantifying the general trend in the overall 
decrease of the sequence became easier. These observations then led us to conjecture 

(13) /3(n) := costt (2,/- +L ) e~ K ^, K = 3.6 ± 0.1, 



as a rough approximation to b n , where L is some real constant. This numerical fit then 
greatly helped us find the main estimate (f3]l above, which has the exact value K = \/2n = 
3.54490 and an additional n 1 / 4 factor modulating the exponential. 

3. The Norlund integral representation 

Our approach to the asymptotic estimation of the b n relies on a complex integral rep- 
resentation of finite differences of an analytic function, to be found in Norlund's classic 
treatise l20l §VIII.5] first published in 1924. In computer science, this representation was 
popularized by Knuth |[T4l p. 138], who attributed it to S.O. Rice, so that it also came to be 
known as "Rice's method"; see |9] for a review. 

Lemma 1. Let 4>(s) be holomorphic in the half-plane 5R(s) > uq — g. Then the finite 
differences of the sequence (<p(k)) admit the integral representation 

fct^ V fc / 2m Jc s(s- l)---(s-n) 

where the contour of integration C encircles the integers {n Q , . . . , rt} in a positive direction 
and is contained in 5R(s) > n — ^. 

Proof. The integral on the right of (fl4l is the sum of its residues at s = no, . . . , n, which 
precisely equals the sum on the left. □ 

An immediate consequence is the following representation for the differences of zeta 
values (S n is as in (0): 



(Choose T > n and consider the finite contour (negatively oriented) consisting of the line 
from § — iT to | + iT, followed by the clockwise arc of |s| = y/T 2 + 9/4 that lies to the 



n \, ("!) 



3/2+ioo 



^The notation x = y designates a numerical approximation of x by y to the last decimal digit stated. 
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right of the given line. The contribution of the circular arc is 0(T ™) and thus vanishes as 
T — > +00.) 

Since b n is 6 n plus a correction term (see Equation (f2]i) and <5„ admits the integral 
representation (0, the first step of our analysis is to move the line of integration further to 
the left. It is well known that the Riemann zeta function is of finite order in any right half- 
plane El §5.1], that is, |C(s)| = 0(\s\ A ) uniformly as |s| — > 00, for some A depending 
on the half-plane under consideration. As a consequence, the integral of ( TTBT l remains 
convergent, when taken along any vertical line left of 0, as soon as n is large enough. 
Under these conditions, it is possible to replace the line of integration 5R(s) — | by the 
line 5ft(s) = — |, upon taking into account the residues of a double pole at s = 1 and a 
simple pole at s = 0. We find in this way 

f-ll n_1 1— 1/2+100 7l' 
Sn = (-l)"" 1 ^! + Ro) + [ ' / COO-? 7T-, \ ds, 



2-7TZ J-l/ 2 -ioo s(s-l)---(s 

where, as shown by a routine calculation: 

(-l) n Ro = -\, {-l) n Ri = -n(l - 7 - #„-i). 
The residues thus compensate exactly for the difference between 5 n and b n , so that 



(-1) 



-1/2+ioo 



(16) b n = v ' / ;ds. 

2tt« J-i/2-ioo s(s - 1) • • ■ (s - n) 

4. Saddle-point analysis of zeta differences 

The integrand in ( TTST l is free of singularities on the left-hand side, and is thus a good 
candidate for estimation by the saddle-point (or steepest-descent) method. This evaluation 
is the main topic of this section, and culminates with the derivation of one of the principal 
results of this paper, TheoremQ] 

We make use of the functional equation of the Riemann zeta function under the form 

(17) C(s) = 2T(1 - S )(2^) s ~ 1 sin y C(l - a). 

From this relation, upon performing the change of variables s 1— > — s, we obtain 

d» «.(?)«!+.) , : ! 1 r . (i+ ," + 

which is the starting point of our asymptotic analysis. 

The integral representation ( fT~8b has several noticeable features. First, the integrand has 
no singularity at all in 5R(s) > h and it appears to decay exponentially fast towards ±ioo 
(see Equations (1201 , ( l2Tb . and ( f2lT i below). This means that one can freely choose the 
abscissa c (with c > |) in the representation 

(19) - 4 si " (?) <(1 + •» .(.fi)'-(. 1> + .) 

The very absence of singularities calls for an application of the saddle- point method. 

The factor £(1 + s) remains bounded in modulus by a constant, and is in fact barely 
distinguishable from 1, as 5R(s) increases, since 

(20) C(^)-l + 0(2^ (s) ), R(«)>|- 
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Also, for large |s|, the complex version of Stirling's formula applies: 

(21) r(l + s) = s s e- s V2™ (1 + 0(|sr 1 )) , 5R(s)>0. 

Finally, the sine factor increases exponentially along vertical lines: one has 



. 7TS [ITS 

2i sin — = — exp —i 



(-i^) + (e-^«/ 2 ) , 9f(a) > 0, 



(22) _ .... 

with a conjugate approximation holding for 3(s) < 0. 

In anticipation of applying saddle-point methods, the approximations d20l l, ( |2TT >. and d22l 
then suggest the function e"^- 1 as a simplified model of the integrand in the upper half- 
plane, where 

n\T(s) 2 



(23) 



u(s) = -slog(27r) - i— 



log 



r( s - 



We shall demonstrate shortly that the location of the appropriate saddle-points in the com- 
plex plane scale as y/n, which may be confirmed by numerical experiments. Therefore, 
in performing an asymptotic analysis, it is appropriate to perform a change of variable 
s — Xy/n, and expand in descending powers of n, presuming x to be approximately con- 
stant. We find, uniformly for x in any compact region of di(x) > 0, ^s(x) > 0: 



(24) 



u){xy/n) = 

u>'(xy/n) = 
U)"(xy/n) = 



2 log x - 2 - log(27r) - 



1 



1 



logn 



log a; + log(27r) - -x 2 + 0(n" 1/2 ) 



log(27r)-i- 
1 -= + 0(n- 1 ). 



2 log x 



(The symbolic manipulation system Maple is a great help in such computations.) 

From the second line of d24l . an approximate root of uj'(s) is obtained by choosing the 
particular value xq of x that cancels u>' to main asymptotic order: 



The substitution of this value 



(25) x = e i7r/4 V27r. 
This corresponds to the following value for s, 

(26) a = cr(n) = XQ\fn = (1 + i)i 

which thus is also an approximate saddle-point for e w< - s '. 
given the first line of (124-b then leads to 

(27) exp (po(a(n))) = exp (2i\/nnj ■ exp (— 2\/nnj ■ 

where IT is an unspecified factor of at most polynomial growth. By using a suitable contour 
that passes though <j(n), we thus expect the quantity in d27l ) to be an approximation (up 
to polynomial factors again) of b n . This back-of-the-envelope calculation does predict 
the exponential decay of b n as exp (— 2y / 7m), in a way consistent with numerical data, 
while the fluctuations, sin + 0(1)) , are seen to be in stunning agreement with the 

empirically obtained formula $13[ . 

We must now fix the contour of integration and provide final approximations. The 
contour adopted (Figure Q]) goes through the saddle- point a = <j(n) and symmetrically 
through its complex conjugate a = a(n). In the upper half-plane, it traverses cr(n) along 
a line of steepest descent whose direction, as determined from the argument of uj" (a), is 
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Figure 1 . Left: The saddle-point contour used for estimating b n . The 
arrows point at the directions of steepest descent from the saddle-points. 
Right: the landscape of the logarithm of the modulus of the integrand in 
the representation of b n for n = 10. 



at an angle of ^ with the horizontal axis. The contour also includes parts of two vertical 
lines of respective abscissae 5ft(s) = c\ \fn and C2\/n, where 

< Ci < y/TT < C 2 < 2y/TT. 

The choice of the abscissae, c\ and C2, is not critical (it is even possible to adapt the 
analysis to c\ = ci = t/tt). One verifies easily, from crude approximations, that the 
contributions arising from the vertical parts of the contour are 0(e _iov/ "), for some Lo > 
2y / 7r, i.e., they are exponentially small in the scale of the problem: 

(28) f =0 (e~ L ^\ L > 

J vertical 

The slanted part of the contour is such that all the estimates of d24l > apply. The scale of 
the problem is dictated by the value of w"(er), which is of order 0(n~ 1 / 2 ). This indicates 
that the "second order" scaling to be adopted is n 1 / 4 . Accordingly, we set 

(29) s = (1 + OVTrn + e 5l7r/8 yn 1/4 . 

Define the central region of the slanted part of the contour by the condition that \y\ < 
log 2 n. Upon slightly varying the value of x around xq, one verifies from ( |24] | that, for 
large n, the quantity 

(xo^ + ^'h^n^j 

is an upward concave function of t near t = 0. There results, in the complement of the 
central part, \y\ > log 2 n, the approximation 



exp 



(lu (xoVn + e 5i7r / 8 yn 1/4 ) ) < e^ x °^ ■ exp (-L x log 2 n) 



Li > 0. 
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Figuratively: 

(30) f = f +0(exp(-L 1 log 2 n)) . 

J slanted J central 

Thus, from ((28) and (13 Ot . only the central part of the slanted region matters asymptotically. 
This applies to but also to the full integrand of the representation ( fT9T > of b n , given 
the approximations (|20l>— (124-b . 

We are finally ready to reap the crop. Take the integral representation of ( fT9b with the 
contour deformed as indicated in Figure Q] and let 6+ be the contribution arising from the 
upper half -plane, to the effect that 

(3D b n = m(b+), 

by conjugacy. In the central region, 

s^x^ + e^^yn 1 ' 4 , 

the integrand of $1% becomes 
(32) 

-i.^) . (27T)- 1 ■ f-lr) ■ fl + 0(2-^)) ■ ^= ■ e ^ xa ^ ] ■ e^ 1 ^ ( 1 + ( -L 



2i J \ ) y/n \ \V™ / 

The various factors found there (compare ( fT9] > to e w ( s ) with u>(s) defined in ( 1231 ) are in 
sequence: the Cauchy integral prefactor; the correction (27r) _1 to the functional equation 
of Riemann zeta; the factor —l/(2i) relating the sine to its exponential approximation; 
the approximation of Riemann zeta; the correction s/(s + n) of the Gamma factors; the 
main term e"'"'; the anticipated local Gaussian approximation; the errors resulting from 
approximations (|20T> — (124-b. which are of relative order 0(n -1 / 2 ). Upon completing the 
tails of the integral and neglecting exponentially small corrections, we get 

(33) b+ = K e"(*°^ i| J + °° e -y 2 /^ dy ■ (W 8 n x / 4 ) (l + O (-^ 

where Ko = — 1/ (47r 2 ) is the constant factor of d32l , while the factor following the integral 
translates the change of variables: ds = e 5 " 1 : ^ 8 n 1 ^ 4 dy. 

The asymptotic form of b n is now completely determined by OTb and d33l . We have 
obtained: 

Theorem 1. The Newton coefficient b n of((s) — l/(s — 1) defined in (f2| satisfies 

(34) b n =[ — \ e- 2 ^cos 

The agreement between asymptotic and exact values is quite good, even for small values 
of n (Figure |2). 

In summary, the foregoing developments justify the validity of applying the saddle-point 
formula to the Norlund Rice integral representation (fT9b of zeta value differences. Under 
its general form, this formula (which, without further assumptions, remains a heuristic) 
reads 



(35) / e-^W dx = a J* - e -W("°) ( l + Q ( 1 



_ Nf"(x ) V \N 

Here the analytic function f(x) must have a (simple) saddle-point at xq, that is, f'(xo) = 
and f"(xo) is the second derivative of / at the saddle-point. In the case of differences of 
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FIGURE 2. A comparative plot of b n and the main term of its approxi- 
mation (|34| |. both multiplied by e 2 ^™n _1//4 , for n = 5 . . 500. 



zeta values, the appropriate scaling parameter is s — x^/n corresponding to N = yjn, and 
the function / is 

f(x) = lim —i=uj (x\fn) , 

up to smaller order corrections that we could treat as constants in the range of the saddle- 
point. As we shall see in Section [6] this paradigm adapts to sums involving Dirichlet 
i-functions. 



5. Convergence of the Newton series of zeta 

The fact that the coefficients b n decay to zero faster than any polynomial in 1 jn implies 
that the Newton series 

oo 

(36) *(s) = £;(-l)% 

ri=0 

with b n given by ©, converges throughout the complex plane, and consequently defines 
an entire function. Set Z(s) := ((s) — l/(s — 1) with Z(l) = 7. We have, by construction 
$(s) = Z(s) at s = 0, 1, 2 ... , but the relation between $ and Z at other points is still 
unclear. 

Corollary 1. The Newton series of (1361 > is a convergent representation of the function 
C(s) — l/(s — 1) valid at all points s G C. 

Proof. Here is our favorite proof. A classic theorem of Carlson (for a discussion and a 
proof, see, e.g., Hardy's Lectures lfl2l pp. 188-191] or Titchmarsh's treatise [23, §5.81]) 
says the following: Assume that (i) g(s) is analytic and such that 

\g(s)\<Ce A ^, 
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where A < tt, in the right half- plane of complex values of s, and (ii) g(0) = g(l) = . . . = 
0. Then g(s) vanishes identically. 

To complete the proof, it suffices to apply Carlson's theorem to the difference g(s) = 
$(s + 2) — Z(s + 2). Condition (ii) is satisfied by construction of the Newton series. 
Condition (i) results from the fact that Z(s + 2) is 0(1) while a general bound due to 
Norlund (Equation (58) of [20, p. 228]) and valid for all convergent Newton series asserts 
that |$(s + 2)| is of growth at most e^N, throughout 3?(s) > -|. □ 

An alternative proof can be given starting from a contour integral representation for 
the remainder of a general Newton series given ll20l p. 223]. Yet another proof derives 
from a turnkey theorem of Norlund, quoted in lfl9l p. 311]: In order that a function F(x) 
should admit a Newton series development, it is necessary and sufficient that F(x) should 
be holomorphic in a certain half-plane ^R(x) > a and should there satisfy the inequality 
\F(x)\ < £721*1, where C is a fixed positive number. In a short note, Baez-Duarte [2] justi- 
fied a similar looking Newton series representation of the zeta function due to Maslanka — 
however his bounds on the Newton coefficients are less precise than ours and his arguments 
(based on a doubly indexed sequence of polynomials) seem to be somewhat problem- 
specific. 



6. DlRICHLET L-FUNCTIONS 

The methods employed to deal with differences of zeta values have a more general 
scope, and we may reasonably expect them to be applicable to other kinds of Dirichlet 
series. Such is indeed the case for any Dirichlet L-function, 



n" 

n=l 



where x is a multiplicative character of some period k, that is, for all integers m, n, one has: 
x(n+k) = x(n), x{ mn ) = x( m )x( n )> x(l) = 1' an d x( n ) = whenever gcd(n, k) ^ 1. 
Let £(s, q) be the Hurwitz zeta function defined by 



oo 



n= + 

Any Dirichlet L-function may be represented as a combination of Hurwitz zeta functions, 

(38) L( X ,s) = 1 £ x ( m ) C ( s , ^ 

where k is the period of \. In particular, the coefficients of the Newton series for L(x, s) 
are simple linear combinations of the quantities 



(39) ^ v ..,.„,-^. t 
which we adopt as our fundamental object of study. 

Theorem 2. The differences of Hurwitz zeta values, A n (m, k) defined by (1391 . satisfy the 
estimate 

(40) A n (m,k) = - - ^ ^(j) +ln* + l--Hn-i +a n (m,k) 
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where the a n (rn, k) are exponentially small: 

1/4 



1 / In \ I Ann \ I Ann 5n 2nm 



a n (m,k) = — I — J exp , 

(41) k \nk / \ \ k I \\k 

+ o(n- 1 / i e- 2 V™^y 
Here ip( x ) = r'(x)/r(x) M f/ze logarithmic derivative of the Gamma function. 

The previous results for Riemann zeta may be regained by setting m = k = 1, so that 

<5„ = A n (l, 1) and 6„ = a n (l, 1). 

Proof. Converting the sum to the Norlund-Rice integral, and extending the contour to in- 
finity, like before, one obtains 

(42) An (m,k)= [ —^-nl / ' ' k > .ds 

2m J|_ j;oo k s s{s - 1) • • • (s - n) 

Moving the contour to the left, one encounters a single pole at s = and a double pole at 
s = 1. The residue of the pole at s = is 



_ x u ( Tn\ 1 m 

Res( S = 0)=c(0, T ) = -- T . 



(See ll26l p. 271] for this evaluation.) The double pole at s = 1 evaluates to 



il>[-r) + \nk+l- H n -i 



Ti 

Res( S = l) = - Lruy 
Combining these, one obtains ( l40b where the a„ are given by 

f— l) n f~i+ io ° C ( s —) 
(43) a „( m , fc ) = i^ n! / ' feZ da . 

As before, the a n (m, k) have the remarkable property of being exponentially small; 
that is, a n (m, k) — O (e~ K ^™^J, for a constant K that only depends on k. The precise 
behavior of the exponentially small term may be obtained by using a saddle-point analysis 
parallel to the one given in the previous sections. Its application here is abbreviated, as 
there are no major differences in the course of the derivations. 

The term a n (m, k) is represented by the integral of (143V At this point, the functional 
equation for the Hurwitz zeta may be applied. This equation is 



(44) 



/ m\ 2T(s) m-^ ( ns 2npm\ ( p\ 

c ( 1 - s '¥J = (2^g cos (T ST J C ( s 'fcJ 



and it either follows from adapting the "second proof" of Riemann for the common zeta 
function 11241 §2.4] or from the transformation formula of Lerch's transcendent <fi found 
in 11261 p. 280]. This allows the integral in d43l to be expressed as a sum: 

a (m k) -™ V r +l °° 1 F(5)r(5 = 1} coc (™ - ^) C (s ds 
a n (m,k)- kriLsJ^ {2n) s T{s + n) ™[ 2 k )<>{ S > k ) ds 
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It proves convenient to pull the phase factor out of the cosine part and write the integral as 

nl ^ f 2irpm\ 
a n (m,k) = — - — > exp i- 

km — ' v 

p=i 



r _J_£(£)E(i^) exp c L 2.) ds + e c 



where c.c. ("complex conjugate") means that i should be replaced by —i in the two exp 
parts. 

To recast the equation above into the form needed for the saddle-point method, an as- 
ymptotic expansion of the integrands needs to be made for large n. As before, the appro- 
priate scaling parameter is x = s/ y/n. The asymptotic expansion is then performed by 
holding x constant, and taking n large. Thus, one writes 
(45) 



a n (m,k) = -t—J2 e i27rpm/fc / e w( - x ^dx +e -^P m / k / 



0=1 

Proceeding, one finds 



(JQ — 2CJO >^<7o — loo 



k + p 



w(s) = logn! + - logn - s log f^l - i — + log r(g + n) + O 

where the approximation \og( (s,p/k) = (k/p) s + O ((p/ (k + p)) s ), for large 3?(s), has 
been made. Expanding to 0(1/ y/n) and collecting terms, one obtains 



ui(xy/n) = — logn — x^/n 
(46) 2 



2tTP 7T 

log— +z- + 2-21oga; 

2 

log 2tt - 2 log X - — + O (n~ 1/2 



The saddle-point is obtained by solving u)'(xy/n) = 0. To lowest order, one has xq = 
(1 + i)y/Ttp/k. Also, ijj"(xy/n) = 2/x^/n + Substituting in the saddle-point 

formula d35l l. one directly finds 
(47) 



e 



while the integral for ZJ is the complex conjugate quantity (having a saddle-point at the 
complex conjugate location). Inserting this into equation (05]) gives a sum of contribu- 
tions for p = 1, . . . , k, of which, for large n, only the p = 1 term is seen to contribute 
significantly. So, one has the estimation fiTI ) of the statement. □ 



7. Perspective 



The previous methods serve to unify and make precise estimates carried out in the litera- 
ture by a diversity of approaches. For instance, the study of quantities arising in connection 
with Li's criterion calls for estimating, in the notations of d39l l. 

(48) A„(l, 2) = £ (") (-l)'(l - 2 "')CW- 
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Bombieri and Lagarias encountered this quantity in |@] Th. 2] and Coffey (see his Si (n) 
in J3) proved, by means of series rearrangements akin to ( [Tol l used in conjunction with 
Euler-Maclaurin summation the inequality 

71 Tl ! 

(49) 4,(1,2) > _logrc+( 7 -l)- + -. 

Our analysis quantifies A n (l, 2) to be 

n 11 
4,(1, 2) = + n( 7 - - + - log 2) + o(l), 

where the o(l) error term above is a n (l, 2), which is exponentially small and oscillating: 

(50) a n (l, 2) = i (-) V4 exp (- V2^) cos f V2^ - ^ ] + (n" 1 /^"^ 



2 Vtt/ "V / V 8 , 

Another observation is that the combination of Norlund-Rice integrals and saddle-point 
estimates applies to many "desingularized" versions of the Riemann zeta function, like 

(1-2 1 - S )C( S ), (s-l)C(s), C(2«) - 2^T> ( 2s - 1 X(2s). 

The first one is directly amenable to Theorem|2] The Newton series involving C(2s) include 
Maslanka's expansion [16| (relative to (2s — l)£(2s)) and have a striking feature — their 
Newton coefficients are polynomials in it with rational coefficients. In addition, the expo- 
nential smallness of error terms in asymptotic expansions of finite differences of this sort 
has the peculiar feature of inducing near-identities that relate rational combinations of zeta 
values and Euler's constant. For instance, defining the following elementary variant of b n , 

we find to more than 35 digits of accuracy, 

c 4 99 -ff 4 99 + l = 0.5782156649 01532 86060 65120 90082 40243... 

7 = 0.5772156649 01532 86060 65120 90082 40243..., 

where the sole discrepancy observed is in the third decimal digit. 

The Norlund integrals are also of interest in the context of differences of inverse zeta 
values, for which curious relations with the Riemann hypothesis have been noticed by Fla- 
jolet and Vallee [ 10], and independently by Baez-Duarte in a scholarly note O. Consider 
the typical quantity 



which arises as coefficient in the Newton series representation of 1/C(s). Its asymptotic 
analysis can be approached by means of a Norlund-Rice representation as noted in related 
contexts by the authors of (10| and more recently by Maslanka in 111811 . The following 
developments provide a rigorous basis for some of the observations made in 1181 . simpli- 
fies the criterion for RH that is implicit in 1101 . and offers an alternative (based on real 
extrapolation and the Mellin transform) to Baez-Duarte's treatment, while pointing in the 
direction of easy generalizations relative to l/£(2s), <^'(s)/((s), £(s — l)/((s), or other 
similar functions. 

Theorem 3. The differences of inverse zeta values d n defined by ( 15 It are such that the 
following two assertions are equivalent: 
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FVBD Hypothesis (akin to EUTO)). For any e > 0, there exists a constant C e > 
such that 

\d n \ < C t n l l 2+ \ 

RH (Riemann hypothesis). The Riemann zeta function £(s) is free of zeros in the 
half-plane 5R(s) > 

Proof, (i) Assume RH. Under RH, it is known that, given any gq > h and any e > 0, one 
has 

(52) -L- = 0(\t\ e ), for R(s) = cr , where t = 9f(s) 

C(s) 

(see Equation (14.2.6) of (24] p. 337]). Then, start from the Norlund integral representation 
(cf Lemma[TJand Equation (15[), 



(53) d n = J„(c), where J„(c) 



(-1)' 



n— 1 /-c+Joo 



27TJ 



1 



C(s) s(a - 1) • • • (a - n) 



ds , 



which is valid unconditionally for c S (1,2). Next, we propose to move the line of inte- 
gration to c = cr - To this effect, observe that the integral J n (ao) defined in d53l converges 
and is O(n <T0 ), since 



|J»(*o)| < 
< 



1 fn-(T 
2-koq \ -cr 

1 fn-a Q 
2nao \ -(T 
0(n CTo ). 



-cr • ■ • (-(T + n) 



(—o"o — it) ■ ■ ■ ( (To — it + n) 
-oo(-0"o + 1) 



(-(T - it)(-£7Q - it + 1) 



There, the second line results from the fact that, for x, t real, one has \ {x/[x — it)\ < l;the 
third line summarizes the asymptotic estimate ("7^°) = 0(n CT ° ) (by Stirling's formula) 
as well as the fact that the integral factor is convergent (since the integrand decays at least 
as fast as 0(\t\- 2+e ) as \t\ -> +oo). 

(ii) Assume FVBD. First, a reorganization similar to the one leading to ( TTOb but based 
on the expansion of 1/ £(s) shows that 



with fi(£) the Mobius function. The general term of the sum decreases like n/£ 2 , which 
ensures absolute convergence. Next, introduce the function 



I>(a0 = 5>(*) 



,-x/i _ 



1 + 



£>1 



whose general term decreases like x 2 /£ 2 . 
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Fix any small S > (<5 = jk is suitable) and define £ — [x 1 s \. The difference 
d n — D(n) satisfies 



d n -D(n) = £>W 



6=1 



1 



1 



_ e -n/e 



(54) 



E + E Me-»" 
K e<£ i>i J 

(e e- n ^) + E G 



a n/^+n log(l-l/£) 



(n 5 ) , 



by series reorganization, a split of the sum according to I ^ £q, and trivial majorizations. 

Given 1541 . the FVBD Hypothesis implies that D{x) = 0(x 1 / 2+e ), at least when £ is a 
positive integer. To extend this estimate to real values of x, it suffices to note that D{x) is 
differentiable on IR>o, and 



= -E 



i 



is proved to be 0(1) by bounding techniques similar to 1541 . Thus, assuming the FVBD 
Hypothesis, the estimate 



(55) 



D(x) 



0(V 



/2+e 



holds for real values of x. 

Regarding the behaviour of D{x) at 0, the general term of D(x) is asymptotic to 
x 2 /(2£ 2 ), so that D(x) — 0(x 2 ), as x — > + . This, combined with the estimate of D(x) 
at infinity expressed by OBI , implies (under the FVBD Hypothesis, still) that the Mellin 
transform 



(56) 



D*(s) 



D{x)x s ~ dx, 



exists and is an analytic function of s for all s in the strip —2 < 5R(s) < — | — e. On the 
other hand, the usual properties of Mellin transforms (see, e.g., the survey |8|) imply that 

T(s) 



(57) 



(oo \ poo 

EmW 8 ] • y o [e 



„s-l 



da; = 



C(-*)' 



at least for s such that —2 < 3?(s) < — 1, which ensures that the usual series expansion 
of 1/C( — s) is absolutely convergent. The comparison of the analytic character of (l56l > in 
-2 < Sft(s) < — 5 — e (implied by the FVBD Hypothesis) and of the explicit form of d57l i 
shows that the Riemann Hypothesis is a consequence of the FVBD Hypothesis. □ 



Numerically, for comparatively low values of n, it would seem that d n tends slowly but 
steadily to 2 = —1/^(0). For instance, we have e?2o = 1.93, c?5o = 1.987, dioo = 1-996, 
^200 = 1-9991. However, it appears from our previous analysis and a residue calculation 
applied to ( f53l ) that there must be complicated oscillations due to the nontrivial zeta zeros — 
these oscillations in fact eventually dominate, though at a rather late stage, as we now 
explain following iflOl . Indeed, assuming for notational convenience the simplicity of the 
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nontrivial zeros of £(s), one has (unconditionally) 

(58) ^EVn ^r+i r( ~f +2 + °^ 

^ C (p) r(n + l-p) 

where the summation extends to all nontrivial zeros p of £(s) with < 'Si(p) < 1, while 
the starred sum (J^*) means that zeros should be suitably grouped, following the careful 
discussion in §9.8 of Titchmarsh's treatise ll24l p. 219] in relation to a formula of Ramanu- 
jan. A simplified model of the sequence d n then follows from the fact that, for large n, any 
individual term of the sum in ( f58l corresponding to a zeta zero p = a + it is asymptotically 

C } C'(p) 

Such a term involves a logarithmically oscillating component, a slowly growing component 
n a (y/n under RH), as well as a multiplier that is likely to be extremely small numerically, 
since it involves the quantity r(— p) x e ^ I I / 2 . For the first nontrivial zeta zero at p = 
i + i 14.13, the term d59l is very roughly 

(60) l(T 9 V"cos (14.13 log n) , 

and for the next zero, at p = \ + i 21.022, the numerical coefficient drops to about 10~ 14 . 
The corresponding oscillations then have the curious feature of being numerically de- 
tectable only for very large values of n: for instance, in order for the first term given by ( f60b 
to attain the value 1, one needs n w 10 18 , while for the contribution of the second zero, 
one would need n w 10 28 . Since it is known that all zeros of £(s) lie on the critical line till 
height To « 5 • 10 s , we can estimate in this way that the presence of a nontrivial zero (if 
any) off the critical line could only be detected in the asymptotic behaviour of d n for values 
of n larger than A*o e 11 " ^q600,ooo : ooo j n summarVi from a numerical point of view, 
a possible failure of RH, though in theory traceable through the asymptotic behavior of d n , 
is in reality violently counterbalanced by the exponential decay of the T factor, hence it 
must remain totally undetectable in practice — analogous facts were observed in iTTOl ITSl . 
It is finally of interest to note that such phenomena do occur in nature, specifically, in the 
determination by Flajolet and Vallee in [ 1 1 of the expected number of continued fraction 
digits that are necessary to sort n real numbers drawn uniformly at random from the unit 
interval. 
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